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Abstract 

The three-dimensional Abelian Chern-Simons theory coupled to a scalar and a fermionic field of 
arbitrary charge is considered in order to study conformal symmetry breakdown and the effective 
potential stability. We present an improved effective potential computation based on two-loop 
calculations and the renormalization group equation: the later allows us to sum up series of terms in 
the effective potential where the power of the logarithms are one, two and three units smaller than 
the total power of coupling constants (i.e., leading, next-to- leading and next-to-next-to- leading 
logarithms). For the sake of this calculation we determined the beta function of the fermion- 
fermion-scalar-scalar interaction and the anomalous dimension of the scalar field. We shown that 
the improved effective potential provides a much more precise determination of the properties of 
the theory in the broken phase, compared to the standard effective potential obtained directly from 
the loop calculations. This happens because the region of the parameter space where dynamical 
symmetry breaking occurs is drastically reduced by the improvement discussed here. 
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I. INTRODUCTION 

Chern-Simons (CS) theory [l| is an important theoretical framework which has been used 
to study many issues on quantum field theory in three space-time dimensions. Among the 
interesting properties of CS theory are the classical conformal invariance and the fact that 
the gauge field does not receive infinite renormalization, leading to a zero beta function for 
the gauge eoupling constant g These are important aspects for the problem of symmetry 
breaking through radiative corrections [3] , which we want to revisit in this work considering 
a CS theory coupled to matter. 

Our study is motivated by some recent developments concerning the summation of the 
power series in the leading and subleading logarithm terms of the effective potential by means 
of the renormalization group equation (RGE) . The RGE allows one to obtain extra 
information from the usual loop approximation, thus providing more refined information 
concerning quantum properties of the model under scrutiny. An important example where 
the RGE have dramatically improved the information obtained in the loop approximation is 
in the analysis of the effective potential for the Standard Model with conformal invariance: 
from the standard one-loop approximation, the effective action of the model does not seems 



to be stable, but with the more precise ap 
it actually is |4|]. For other examples see 



iproximation obtained using the RGE, one discover 

We show here an improved calculation of the effective potential of the theory of a CS 
field coupled to scalar and fermionic fields. The computation includes infinite summations of 
terms of the effective potential which can be carried out with the RGE and the knowledge of 
the elements figuring in it: the beta functions, scalar field anomalous dimension and the first 
logarithm corrections for the effective potential. These elements at lowest approximation 
need a two-loop calculation to be determined since there are no one-loop divergences in odd 
space-time dimensions when using the sort of regularization we adopt here (regularization 
by dimensional reduction). In fact all one particle irreducible diagrams with an odd number 
of loops will be finite under this scheme. Some of the needed elements were computed in 
Refs 2j, l9Ml2|: in this work, we calculate the fermion anomalous dimension and the beta 
function for the Yukawa coupling. 

A peculiarity of CS theory has to be mentioned at this point. The theory involves the 
Levi-Civita tensor which cannot be easily extended to arbitrary dimensions as needed in the 



2 



context of dimensional regularization. A regularization procedure called dimensional reduo 



tion 



13| has been shown to be appropriated in dealing with CS theory 111] : essentially. 



it consists in performing the tensor and gamma matrices algebra in three dimensions, and 
extending only the momentum integrals to arbitrary dimensions. 

The two-loop results, in conjunction with the RGE, allows us to sum up all terms in 
the effective potential where the total power of the coupling constants is one and two units 
larger than the power of the logarithms log (<f)/fi 2 ) (called leading logarithms, LL, and next- 
to-leading logarithms, NLL, terms), as well as some subseries where they are three units 
larger (the next-to-next-to-leading logarithms terms). We study the dynamical symmetry 
breaking of the conformal symmetry in this theory, showing that the improved effective 
potential leads to a much finer determination of the properties in the broken phase, such 
as mass and coupling constant of the scalar field. This happens because the region of the 
parameter space of the theory, where the dynamical breaking of symmetry is operational 
at the perturbative level, is much smaller when considering the improved effective potential 
than for the initial two-loop potential. Another interesting aspect is that, for certain values 
of the parameters, we found two broken vacua, which leads to different physical properties. 
This happens both for the improved and the original effective potential, but the region of 
the parameter space where this happens is much more restricted for the former case. Again, 
the improvement of the perturbative effective potential calculation provides more precise 
determination of the properties of the theory. 

We believe that the outcomes of our analysis involving the Chern-Simons theory en- 
forces the idea that one has to extract the maximum amount of information from a given 
perturbative calculation, by using the renormalization group equations to obtain a better ap- 
proximation to the effective potential. Even if its natural at a first moment to use one-loop 
results to predict masses and coupling constants from any of the many proposed exten- 
sions to the Standard Model, for example, one should enrich the analysis of the dynamical 
symmetry breaking by means of the RGE. 

This paper is organized as follows. The method of using the RGE to sum up series 
of perturbative corrections to the effective potential is outlined in Sec. [Til The model we 
shall study is described in Sec. lIIII Technical details of the two- loop calculations needed for 
this work are presented in Sec.UVI Sec. [V] contains the detailed calculation of the improved 
effective potential, which is used to study the dynamical breaking of the conformal symmetry 
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in Sec. lVIl Finally, our conclusions are summarized in Sec. lVIII 



II. GENERAL CONSIDERATIONS 



We start by reviewing the use of the RGE to calculate the improved effective potential. As 



discussed in 



15] 



the standard practice for solving the RGE by replacing the couplings in the 
effective potential by their running values amounts to a particular application of the method 
of characteristics to solve partial differential equations. This procedure does not exhaust, 
however, the information that is contained in the RGE: actually, a finer approximation can 
be obtaining by writing the effective potential as a general power series in the couplings and 
logarithms of the scalar field, and using the RGE to sum up some infinite subsets of this 
power series. 

To explain the general procedure, we will consider a general model of a scalar field 99 
with a self-interaction of the form ip N , together with interactions with other dynamical 
fields. As known, in three spacetime dimensions, renormalizability imposes that N < 6, but 
we shall not fix any particular value of N in this Section. Let A = {Aj, i = 1, . . . M} denote 
collectively the set of all coupling constants of the theory. The RGE for the regularized 
effective potential V q q (0) reads 



d n d . d 



V eff (0;/i,£,A,L) = O (1) 



(in this section, the sum over all Aj will always be implicit). Here, \x is the arbitrary mass 
scale introduced when we use dimensional regularization to extended the theory to dimension 
D, 7^ is the anomalous dimension of the scalar ip, e = 3 — D, 

L = ln^, (2) 

and 4> is the vacuum expectation value of the scalar field (p. 

For the sake of convenience, we introduce the notation (from now on, we omit the explicit 
dependence on the parameters /i, e, A, L), 

V eS (0) = 4> N S eE (0) , (3) 

where S" e g(0), on very general grounds, is a sum of terms involving different powers of A 
and L, which in principle can be calculated order by order in the loop expansion. 



In order to use the RGE we shall organize the terms in £ g (0) according to the power 
of L relative to the aggregate powers of the couplings A, i.e., 



NLL 



^N2LL 



where 



n>l 

is the sum of the leading logarithms in S q q (0), and 

^e N ff L ^) = E^ NLLAnLn_: 



n>3 



5 eff LL (<P) = J2 C n 2LLxnLn ~ ; 



(4) 



(5) 



(6) 



(7) 



n>3 



are the next-to-leading and next-to-next-to-leading logarithms terms, respectively; here, A" = 
Y[ A"" with rij = n. The RGE allows one to calculate these sums once their first coefficient 
is known, if we have enough information on the /3-functions and the anomalous dimension 
of the scalar field. To see how this come about, we use the definition ([3]) in Eq. flTJ, and take 
Eq. (T5]) into account to rewrite the RGE in a more convenient form, 



S eS ((f); fi, e, A, L) = 



We shall write 7^ and 3\ in the form 



X 2 ) 



„(3) 



+ 



7^+7^ 
/3f + 0? + ■ 



(8) 

(9) 
(10) 



where jjp and 8^ denotes the terms of order A J of the anomalous dimension and beta 
function, respectively; these can be obtained by explicit loop calculations. 

Substituting the expansion (TJJ in (jSJ) we find, at the leading order (terms proportional to 



d 



dL + dX 



0. 



(11) 



This results in a first order difference equation for the coefficients C^ L ; in this way S^((j)) 
can be determined once we know 8^ and the initial coefficient CY". Having S^l. ((f)) at our 



disposal, we can focus at terms of order X n L n 3 in ([8 



-— + {2) — 



5 e L ff^ 



0. 



(12) 
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Since S^(<p) is known, this equation allows us to calculate £gg L (0) if we have (3 X , 7^ 
and C 2 NLL . 

This procedure can be repeated until we have exhausted the information on j3\, 7^ and the 
initial coefficients C from the explicit loop calculations. In summary, the RGE allows one to 
use the knowledge of S^g (0), (3\ and 7^ up to a given loop order to sum up complete subsets 
of contributions for the effective potential arising from all loop orders, thus extracting the 
maximum amount of information from our perturbative calculation. 



III. THE MODEL 



We shall now consider a Chern-Simons field A p in three spacetime dimensions coupled to 
a two component Dirac field ip and a complex scalar field ip, both charged under the U (1) 
gauge symmetry of the CS field according to the Lagrangian 

C =^ p A»d»AP + ifrfD^ + (D»<p? (Dtf) (13) 
~ g (v?V) 3 ~ ay^W^- 

The theory has a self-interaction for the scalar field and an Yukawa-like interaction between 
scalar and fermions fields. In Eq. ( fl"3l) . v is a positive coupling constant and _D M = d^—ienA 11 , 
where n is the charge of the field _D M is acting on. Without loss of generality, we can consider 
n 9 — 1, since any n 9 ^ 1 can be reabsorbed by a redefinition of the gauge coupling constant 
e. Therefore, we will denote simply by n the charge of the fermion, from now on. The 
spacetime metric is g^ u = (1,-1,-1), the fully antisymmetric Levi-Civita tensor e pup is 
normalized as eoi2 = 1, and the gamma matrices were chosen as 7^ = (a 3 , ia 1 , ia 2 ). 

The Lagrangian in Eq. ( fl3l) is a (2 + 1) dimensional analog of the well known Coleman- 
Weinberg model in (3 + 1) dimensions [3j, in the sense that all parameters appearing in the 
classical Lagrangian are dimensionless, so it posseses classical conformal invariance. As we 
assume such an invariance at the classical level, to deal with quantum corrections it is appro- 
priate to use a regularization method that violates it minimally [l^] . The observations made 
in [jjj] regarding dimensional regularization are straightforwardly generalized for regulariza- 
tion by dimensional reduction, which has been used to obtain the quantities we need here. 
Divergent integrals are regulated by the replacement J d 3 k/(2ir) 3 — > // j d 3 ~ e k/(2ir) 3 ^ e , 
where the mass scale /x is introduced to keep the dimensions of the relevant quantities un- 



changed. Conformal invariance is broken explicitly by this mass scale, but /i comes with 
the evanescent exponent e and this, in conjunction with the poles 1/e, means that \i always 
appears inside a logarithm. Also, regularization by dimensional reduction has been shown 
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15|- 



to preserve Ward identities at least until the two loop order 

Details of the two-loop calculation of the effective potential for a theory like in eq. ( I13p 
can be found in [12[. In summary, after introducing a convenient gauge fixing, one defines a 
Lagrangian £j n ^ shifting the scalar fields by a constant, and disregarding terms independent 



of or linear on the fields 
of 



161 ] : after that, the effective potential can be calculated by means 



v , 6 i 



d 3 k 



In [det (iA~l (Jfe,&))] + « < 0|Texpi / d 3 xC int \0 > . (14) 



48" 2 J (2tt) 3 

Hereafter, 2 stands for y/2 (ip^<p}. The first and second terms in Eq. (114p are, respectively, 
the tree approximation and the one-loop correction to the effective potential; the third term 
is the sum of the vacuum diagrams with two and more loops. 



it 



We quote here the two-loop effective potential in the following form 
where (0) is more conveniently written in terms of the coupling constants 



(15) 



v 
8^2 



X 



2tt 



a 
2^ 



(16) 



as follows, 



% (0 = y+ 



n: 



33 



x 



x 2 y + 14y 2 + 



24 ( I 

15 2 ., 1 /• 2 :; > 2 2 

— yz — 6z — bn^x z + 3n^x z 



L. 



(17) 



On the other hand, as discussed in Section|TTl the general form for £? g (0) can be cast as in 
Eq. OH), with 



S lk^= E C^x-y^L^+o- 1 

n,m,p>0 
(n+m+p>l) 

n,m,p>0 
(n+m+p>2) 



(18) 



(19) 



and 



SS L (<P)= £ C^y y m z?L n+m - 3 . (20) 



n,m,p>0 
(n+m+p>3) 



It is known that the beta function of the gauge coupling x vanishes in CS model coupled 



to scalar and fermionic fields j2|; we calculate the two- loop approximation the beta function 
(3 a of the Yukawa coupling, as well as the scalar anomalous dimension 7^ in SectionHVl The 
Renormalization Group equation reads, in our model, 



SeffW>)=0. (21) 



By following the procedure outlined in Section [TTJ, we obtained closed-form expressions for 
S'eff (0), £gjj L (0) and 5g| LL (0). The technical details of this calculation are quite involved 
and are developed in Section|V] The results we obtain are the following, 

S%{*) = %, (22a) 
Sg L (0) = (w) + z 2 S^ (w) , (22b) 

S N e £ L (0) = (x 4 S?^ (w) + z*S^ (w) 

Wz'Sf 2 ^ (w) + x 3 zS^ (w)) L , (22c) 

where 

w = 1 — a x yL (23) 
and the functions of w appearing in Eq. (1221) are explicitly displayed in SectionlVl 



IV. TWO-LOOP WAVEFUNCTION RENORMALIZATION AND FUNCTIONS 

For the purposes of this work we need to calculate the beta function for the Yukawa 
coupling a<p^ \pipijj , which implies in calculating the renormalization of the four-point (p'tpipx/j 
function, as well as the wave function renormalization of the ip field. To evaluate these 
quantities, we calculated in the two-loop approximation the divergent parts of the fermion 
two-point vertex-function IV, and the four point vertex function F t 7^. Free propagators 
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for fermionic, scalar and gauge fields are given respectively by 

A * {k) = ' (24a) 

A„W = j^-, (24b) 
V(*) = ^. (240 

while the elementary vertices are 

trilinear i/ji/jA p -h- —ie^ p ^ £ ^ 2 , (25a) 

trilinear (p) <£> (—q) A p -H- — ze (p + g) //' 2 , (25b) 

quadrilinear Tp<pA p A a -H- ie 2 g pa [i e , (25c) 

quadrilinear ij)tj)Tp(f <-> —ia\i £ , (25d) 

where, in the vertex, the indicated momenta are the ones entering the respective line. 

The diagrams involved in calculating the two-point vertex function of the fermion are 
shown in Fig.HJ and the corresponding divergent parts are given by, 

(«) = -y, (&) = (e) = |^, (26a) 

(<0 = -Y> ( e ) = y. ( 26b ) 

apart from an (ifi 2e /16tt 2 ) ft/e factor. 

We also evaluated the divergent part of the four-point ipipTpip vertex function in the 
two-loop approximation. Our method for this calculation was the following one: all two- 
loop 1PI diag rams for such vertex function were generated using the Mathematica package 



FeynArts 13], resulting in about 200 diagrams The identification of the divergent diagrams 
was greatly facilitated by the fact that, for the purpose of evaluating the divergent part of 
the i/jipTptp function, we could calculate the diagrams with vanishing external momenta. This 
allowed us to prove an important rule, all diagrams with a trilinear TptpA vertex attached 
to an external line are finite due to the antisymmetry of the gauge propagator A^„. This 
rule is graphically represented in Fig.HJ There are also some one-loop diagrams that vanish 
(those depicted in Fig. E]) and appear as subdiagrams of some of the initial set. Using the 
pattern-matching capabilities of Mathematica, we could use such rules to narrow down the 
set of possibly divergent two-loop diagrams to those appearing in Fig. HJ The result of the 
calculation of these diagrams appears in Table [D 



With these results, we can now write down the relation between bare (denoted by the 
subscript zero) and renormalized fields and coupling constants 



y?o = Z*<p = (1 + A)*<p, 

^ = zl^ = (l + F)h, 
a Q Tp Lp Q il) ^Q = ii £ (a + 8a) Jpipipip . 



(27a) 
(27b) 
(27c) 



The constant Z v has already been calculated in 12] , and the results of Eq. ( )26l) and TableU 
allow us to find and 8 a: 

1 



8a 



32tt% 



(7e 4 a + 4eV + 20e 6 - 4a 3 ) 



1 + 



1 



288tt% 



1 



16n 2 e 



(3a 2 - 8e 8 ) 



3-P + *.")-„ 



(28a) 
(28b) 
(28c) 



The beta function for the Yukawa coupling is calculated from the relation (|2Tb ) 

jjf (a + 8a) 



G 



Z^Z^ 



(29) 



since dao/dfi = 0, we have 

da 



dfj, 8tt 2 



Of " -[ ( ^ + ^- ) ae A + a 2 e 2 - a 3 
3o 



(30) 



In terms of the rescaled coupling constants in Eq. ( TT6|) . 



Pa 



97 n 2 \ 2 1 2 



, 21 H — z x z 

72 6 7 2 2 



(31) 



2tt 2 

From Eqs. (128b) and (128b). we obtain the anomalous dimensions for scalar and fermion 
fields, 

1 dZ t , 



7 V 
7</> 



7 n 2 \ 2 , 1 _ 2 



1 dZ, 



Z^p fj, 



1 x H z 

12 6 / 24 



1 2 * 2 

— ar H z 

9 24 



(32a) 
(32b) 



where 7^ has been quoted just for completeness. 

As for the beta function of the coupling u, it is most easily calculated by relating it 
with the effective potential in Eq. (JT7|) and the anomalous dimension 7^ by means of the 
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renormalization group equation, as done in [12|. Here, we just quote the result, taking into 
account Eq. ffT6l) and the fact that the fermion has charge n, 



P v =24 ^1 + ^ x 4 - (n 2 + 20) x 2 y + Uy 2 
+ Ayz 2 - 3z 4 - Qn 2 x 3 z + 3n 2 x 2 z 2 . 



(33) 
(34) 




y 



(a) 



(r/l 



(/) 



Figure 1: Two- loop contributions to the fermion two-point vertex function. 




= finite 



Figure 2: A simple rule for establishing the finiteness of a subset of diagrams: since the external 
momenta can be taken to zero, whenever there is a trilinear TpipA vertex attached to an external 
line, the resulting Feynman integrand would contain a factor x (—iep^), thus vanishing due 
to the antisymmetry of the gauge propagator A^ u . 
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Figure 3: One-loop vanishing diagrams that appear as subgraphs of some of the two-loop contri- 
butions to the four-point vertex function. 
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D6 D7 D8 D9 DIO 




D21 D22 D23 D24 D25 

Figure 4: Potentially divergent two-loop diagrams. 



Dl 


-fae 4 


D6 


-e 6 


Dll 


±a* 


D16 


4 

— ae 


D21 


-e 6 


D2 





D7 


4 

ae 


D12 




D17 


|aV 


D22 


— \ae A 


D3 


-a 2 e 2 


D8 


^ae 4 


D13 


— ^ae 4 


D18 





D23 


4 

— ae 


D4 


4 

ae 


D9 


e 6 


D14 


-±a 2 e 2 


D19 


-|ae 4 


D24 


-e 6 


D5 


^ae 4 


DIO 


I« 3 


D15 


-2e 6 


D20 


-|ae 4 


D25 


-e 6 



Table I: Divergent parts of the diagrams appearing in Fig. 01 omitting an overall factor of i/j, 6 /8ir 2 e. 
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V. CALCULATION OF THE IMPROVED EFFECTIVE POTENTIAL 



In this Section, we apply the methodology outlined in Section|TT] to the present theory. 
We use as a starting point the two-loop effective potential in Eq. (|17[) . from which one can 
identify the numerical values of the initial C m ^ p coefficients of the expansion 

n,m,p>0 
n,m,p>0 

+ E C^ L p x n y m z p L n+m+p - 3 + • ■ ■ , (35) 

n,m,p>0 

by casting Eq. (|T7|) as 

Sf s (<P) =V (C$, + CftjyL) + {x 2 C^yL + z 2 C^) yL 

+ (C L - 4 + C L 4 L ^ 4 + C™tJz + C™fx 2 z 2 ) L . (36) 



-d + 2 7 ,)A + ^| + ^_ 67v 



The beta functions and anomalous dimension that appears in the RGE, 

S eS (0) = ; (37) 
were presented in Section HVl and can be cast as 

Pv = W + ^ + 4 4) . (38) 

where 



/3< 2 > = a W 2 ; 0® = a 2 x 2 y + 



a 3 yz 



'y "xa j 

= a^x 4 + a^z 4 + a§x 3 z + a^x 2 z 2 , (39) 

and 

g = (3^ = b lZ 3 + b 2 z 2 x + b 3 zx 2 + hx 3 , (40) 
as for the anomalous dimension, we have 

7^ = 7? } = cix 2 + c 2 z 2 . (41) 
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The numerical values of the coefficients appearing in the last equations are 



33 

a x = 14, a 2 = 6ci - — = - (n 2 + 20) , 

15 " / n 2 \ 

a 3 = 6c 2 + — = 4, a 4 = 24 M + — J , 



Cl 



+ 



— 6n 2 


a 7 = 


: 3n 2 . 


n 2 \ 




1 




c 2 = 


24' 




97 


n 2 










~ 72 


+ T 



(42a) 
(42b) 
(42c) 
(42d) 

(42e) 



where n is the charge of the fermionic field. 

Using these results, we can split Eq. ( |3"7|) according to the relative powers of coupling 
constants and logarithms, 



dL Py dy 



qLL 



+ 
+ 



l y dy ' z dz lip 



S ek + 



-— + ^)- 
dL v dy 



oNLL 

^eff 



lx dL f y dy 



f y dy ' z dz 7<fi 



C NLL | 

^eff + 



d_ 
dL 



d 



■— + /3 (2) — 
" Py dy 







C N2LL 

^eff 



(43) 



A. Leading logarithms 

Focusing first on terms of order x m y n z p ]_, n + m +p- 2 [ n Eq. fj4*3l . one obtains 

-iLL , 

'eff' 



dL > y dy 



(0) 



d 



d 



dL dy 



C^ p x m y n zVL m+n+ v~ l = 



(m+n+p>l) 



which furnishes the following relation for the coefficients 

{m + n + p - 1) - (n - 1) aiC£ n _ ltP = (m + n + p > 2) 

We find convenient to recast S^g (0) as 

Sjfffo) =^(o) («) + E ^i m+P_1 ^) («) . 

m+p>l 



(44) 



(45) 



(46) 
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where 

5 ( L o)W=EOi,o^ (47) 

n>0 

S(m,p) ( u ) = C^ n p u n , (48) 

n>0 

in terms of the variable u = yL. Inspection of Eq. ( 1T7|) allows one to find the initial coefficient 
of these sums, 

riLL _ i (-iLL _ 
u 0,l,0 ~~ °0,2,0 ~~ "1 ) 

CLL s~iLL s^iLL /^iLL /~iLL /~iLL /~iLL r\ fAQ\ 
1,0,0 — °0,0,1 — °1,1,0 ~~ °1,0,1 ~~ °0,1,1 ~~ °2,0,0 ~~ °0,0,2 — u ■ l^ y J 

By looking at Eq. (145 p with m = p = 0, we have 

C = fl iCi,o (n>2), (50) 
with, together with the values Cq f , C^fo from Eq. fT49l . leads to 

Co L n,o = «r 1 (n>l) (51) 

hence 



5 ( L ] )H = E a > n = r^- ( 52 ) 

n>0 1 

Now setting m = 1 and p = in Eq. (1451) , 

"Clio - (n - 1) ai^w,o = (n > 1) , (53) 

and from this equation one concludes that C\\ Q = 0, which is consistent with the results 
obtained from the two-loop calculation of V^g in Eq. ( )49l) ; this is an important consistency 
check of that result. Also from Eq. (|53p . by recurrence we have 

0^ = (n>0), (54) 

so that S^q) (u) = 0. Similar results are found by setting m = and p = 1, i.e., 



QS,i = (n>0), (55) 



thus ^ L 1} (u) = 0. 
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Now looking at the terms with m + p > 2 in Eq. (1451) . for n = 1 we immediately obtain 

Cgk, = , (56) 
which, by recurrence for larger n, implies that 

Ct L n , p = 0. (m + p>2) (57) 

Summarizing this results, 

Sf^ p) (u) = (m,p^0), (58) 

therefore, 

S% (0) = yS% (u) = V - , (59) 



where we have introduced the definition 



w = 1 — a\u = 1 — aiyL . 



(60) 



B. Next-to-leading logarithms 



Having found S^q, we can now consider terms of order x m y n z Pi J m + n +p 3 j n Eq. (T43 
d „ m d 



8L Py dy 



^eff ~ U - 



(61) 



At this point, the first term is completely known, and we proceed to find out $gg L which, 
as before, will be written in the form 



S^ L (0) =y 2 S^ L (u) + x m z P L m+p - 2 S { 

m+p>l 



NLL ( \ 
(m,p) \ u ) i 



D (0) \ U ) ~ 2^1 U m,n+2,p U 
n>0 



(62) 



(63) 



E^SX- (64) 

n>0 

After some manipulations, Eq. (IBTj) can be cast as 
£ ( H2 - 6 Cl ] x 2 + [na 3 - 6c 2 ] ^} C^y^' 1 

J2 [-(m + n + p-2) Cl% + 01 (n - 1) C^lJ x m y W"^" 8 



n>l 



m,n,p 
(n>l,m+n+p>3) 



(65) 
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Some initial coefficients for these sums are obtained from Eq. (1171) . as follows, 

u 2 

/-tN l^Lj /-tiy (~1L\ /-tlVLLj /—tl\ L/L/ /-tiyLL r\ 

Ogn — Lrnn — >- y nn r > — ^nn — t^ini — >- / 9nn — U 
r l 



tNLL 
210 


33 

~ ~T' 


/-tNLL 
°012 _ 


15 

~ T' 






tNLL 
Oil 


/-tNLL 
°020 


r NLL 
— °002 


/-tNLL 
°110 


/-tNLL 
— °101 


/-tNLL 
— °200 


tNLL 
021 


/-tNLL 
— °003 


/-tNLL 
°030 


r NLL 
— Ly 102 


/-tNLL 
— °111 


/-tNLL 
— °120 



/-tNLL 
°201 



/-tNLL 
°300 



0. 



(66) 



As before, we look at some particular subseries in Eq. ( 1631) . First, isolating terms with 
m = 2 and p = 0, 



[na 2 - 6d] C % - nC«^ + a x (n - 1) C 2 N n L \ = , 
which is consistent with the coefficients found in Eq. ( 1631) . since 

[a 2 " 6cx] C L to - = 6d - f - 6 Cl - (-f ) = 



(67) 



(68) 



Also from Eq. ( 1671) . multiplying by u n 1 and summing up over n, we obtain a differential 



equation for the function S 



NLL 

(2,0) 



U), 



j cLL J cNLL J cNLL 

flD (0) . / c \ oLL a °(2,0) aO (2,0) „ 

a2U ~du~ + (a2 - 6Cl) ^ aiU ^T = ° ' 

or, rewritten in terms of the variable w = 1 — aiU, and according to Eq. (1531) . 

^(2^ «2 6ci 



(69) 



CLxW- 



dw w 



w 



(70) 



The solution can be found satisfying the initial condition S'ffik (w = 1) = as 



oNLL _ °2 



— - 1 



6ci 

a>i 



(2,0) 
1 

1 

w 



Proceeding similarly for terms with m = and p = 2 in Eq. ( 165]) . we have 
K - 6c 2 ] - nCgJ + oi (n - 1) C N ^ 2 = , 
whose consistency with the initial values in Eq. (163|) can also be checked, 



[03 - 6c 2 ] C$j - C» LL 



"012 



6c 2 + — - 6c 2 - — = . 



(71) 



(72) 



(73) 



Eq. (172]) furnishes a differential equation for S^J^ whose solution is 

a 3 



oNLL _ -o 

5 (°< 2 ) " 237 



1 

— - 1 



6c 2 
a i 



1 

1 



(74) 
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For all remaining terms in Eq. (165]) . the relation 

-(m + n + p-2) C^ p + en (n - 1) C^_ lp = , (75) 

together with the initial coefficients C^ L L p , and C^ L 2 L p in Eq. (IMjl , implies that C^ L n L p = 

0. 

This way, the only nonvanishing subseries of $gg L (0) are the ones defining Sjf^ and 
^(02)' anc ^ we enc ^ U P with 

5g L (0) = + z 2 ^ . (76) 



C. Next-to-next to leading logarithms 



Finally, we focus on terms proportional to x m y n z p L m+n+p 4 
d „™ 9 



Jx dL Py dy 



/ 5(3)_^ +/3 (3)_^_ 6 (2) 

+ h n ' — 

dL y dy 



oNLL 

^eff 



oN2LL _ n 

^eff - U ' 



(77) 



This time we only have information from the two-loop computation of V^g of the following 
initial coefficients, 

/ r 2 ' 

— . , .,, , — , ,,, , jL — OA I 1 _L _ 

y 004 ~~ °> ^202 ~~ tJlb ) ^301 ~~ u "' ' ^400 



(78) 



so we will focus on the subseries of terms of the form x A y n L n , z 4 y n L n , x 2 z 2 y n L n , and x 3 zy n L n 
in Eq. (EZD. 

We start with terms proportional to x 4 y n L n ; from Eqs. fl38|) to fj4Tl) . Eqs. (1591) and (TTSD . 
they arrive from the following terms of Eq. (I77|) , 



aix 



d_ 

dy 



5 (o) + 



2 5 fi 2 

^2^ ?/7r ocia; 

dy 



2 oNLL , 
X D (2,0) + 











oN2LL _ n 

^eff ~ U 



(79) 



or, writing explicitly, apart from the overall x 4 factor, 

£ W^ 1 + - 6c i) cJS^] + £ [- ( n + x ) + «i^4 N n L o L ^ n+1 ] = • 

n>l n>0 

(80) 

This relation is consistent with the initial coefficients in Eq. ( 179 p . since for the term propor- 



tional to u° we have 



a 4Cai,o — — 04 — 24 I 1 + — 



yN2LL 



ir 



(81) 
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From Eq. ( 18 Op we obtain the relation 



(n + 1) a 4 C^ n+h0 u n + (na 2 - 6c,) C^u n -{n + 1) C™fu n + a 1 nC^u n+1 = 



which provides the following differential equation 



° 4 { u iL + x ) s $ + { a2U ^i ~ 6ci ) s w> + ( M {aiU ~ 1} I: ~ x ) s ™v = 



to be solved for 

oN2LL /A _ V^ r N2LL n 
D (4,0) W — 2-*i 4,n,0 M • 
n=0 

Eq. ( 18 3 p is more easily solved when written in terms of the variable w = 1 — aiyL, 
w (w - 1) 4- + A C - f«2 («; - 1) 4~ - 6 Cl ) Sf LL 



+ a 4 [(w-l) — + l)S$. 



The solution 5gf L } L is 



Q N2LL _ «3 , «2 «1 

y ' w w 



where the coefficients «j are 



a 3 



_a|_ 
4a x 



3a2Ci 

a 2 = 1 , 

a\ Via, 3 

18cf an 04 

ai liai 3 

3a 2 C! 18cf a| 04 

a, d\ Yld\ 3 

Proceeding similarly for terms of the form z A y n L n , we obtain the relation 

(n + 1) d 5 C^ L n+lfi u n + (nd 3 + 2h - 6c 2 ) C^u n 

-(n + l) C^u n + fll nCV +1 = (n > 0) 

which provides us a differential equation for the determination of 

C.N2LL _ /^N2LL„ n 
(0,4) ~ 2.^1 0,n,4 " 5 
n=0 
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as follows, 

d _ \ ^ T t f d 



a 5 ( u— + 1 ) Sfk + ( a 3 u^- + 26i - 6c 2 ^) + (a lM - 1) A _ ^ 5( N2LL = Q (90) 



The solution, again in terms of the variable w, is 



where 



a 



2 



& = ^ , (92a) 
4ai 



^2 = 



03&i _ 3a 3 C2 _ Ji3_ + 05 , g2b x 
3ai ax 12ai 3 ' 
66ic 2 + 0361 + 18c| _ _a|_ + 05 



?o = 



ax 3ai ai 12a! 3 ' 
66ic 2 2a 3 6i + 3a 3 c 2 _ 18c 2 _ aj + 05 



ax 3ai ai ai 12ai 3 
Now, focusing on terms proportional to x 2 z 2 y n L n , we obtain the relation 



(n + 1) a 7 C^ +li0 ^ + (a 2 C N S + a 3 ^ L ) nu n + (2o 3 - 6 Cl ) C^u n 

-6 C2 C 2 N n >" " (n + 1) C™V + ai^2 N SN" +1 = ■ ( 93 ) 

The function 

^ = £^ N £N n , (94) 



n=0 



is determined by the equation 



+ (a 3 «^ - 6c 2 ) Sgft + (am - 1) £ - l) S ( N 2 2 2 L } L = . (95) 



whose solution is 



Q N2LL 73 72 7l , o/ / Qfi x 

^(2,2) = — 3 + — H ^7o, (9o) 



where 



73 = ^ , (97a) 

a 3 o 3 3a 3 ci 3a 2 c 2 a 2 a 3 a 7 

72 = ^— + ' ( 97b ) 

3ai ai ai oai 3 

6o 3 c 2 a 3 o 3 36cic 2 a 2 a 3 a 7 

7i = + ^— + -7, — + t , (97c) 

01 oai ai oai 3 

66 3 c 2 2a 3 o 3 3a 3 ci 3a 2 c 2 36cic 2 a 2 a 3 a 7 

7o = 5 + + ~ a — + -5- • (97d) 

ai oai ai ai ai oai 3 
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Finally, summing up terms of the form x 3 zy n L n , we have the relation 

(n + 1) a 6 C^ n+lfi u n + 2b,C^ 2 u n - (n + 1) C«" + «i n C« n+1 = , (98) 
which determines 

CN2LL _ /^<N2LL 
D (3,l) -Z^°3,n,l 

by the equation 



. NT2LL^n 



n>0 



d 



dw 



a 6 ( + 1 ] S$ + 2b 4 S L ^ + ( u (a lU - 1) - 1 ) 5, 



-iNLL 
'(0,2) 



d 



dw 



~iN2LL 
'(3,1) 



The solution reads 



with 



C N2LL _ ^2 5i 
^(3,1) — ~~ o" H rOo, 



3ai 3 ' 
66 4 c 2 a 3 6 4 , a 6 

Ol = 1" T, 1" ~7T , 



s 



a\ 3ai 3 
664C2 2a 3 6 4 a6 

ai 3ai 3 



As a result, 



C-N2LL _ I 4qN2LL , 4 qN2LL , 2 2rtN2LL , 3 qN2LL\ r 

eff ~~ v ( 4 >°) + (°> 4 ) + ( 2 . 2 ) + 2,:, (3,i) J ^ • 



(99) 



(100) 
(101) 

(102a) 
(102b) 
(102c) 

(103) 



VI. DYNAMICAL BREAKING OF SYMMETRY 



In this section, we show how the dynamical breaking of conformal symmetry occurs in 
the present theory, taking into account the improved effective potential we have obtained, 



^eff (0 =\** K L ff (*) + SeT (0) + S e N ff L (0) + «} , 



(104) 



« being a finite renormalization constant, which is determined by imposing the tree level 
definition of the coupling constant 



dV eff (0) 



gnreeW 
d0 6 



= 6!7r 2 2/ • 



(105) 



The fact that V e g (0) has a minimum at 2 = \i requires that 

dV eS (0) 



d0 



o, 



(106) 
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Figure 5: Sections of the parameter space of constant e 2 , showing where the dynamical symmetry 
breaking occurs, using the improved effective potential. 

and this equation is used to determine the value of y as a function of the free parameters x, 
z and n. This give us a seventh-degree equation in y, and among its solutions we will look 
for those which are real and positive, and correspond to a minimum of the potential, i.e., 

.2 _ dXtti^) 



m, 



> 



(107) 



03 =A « 



We explore the parameter space of the constants x, z, n, looking for values where the 
dynamical symmetry breaking is operational at the perturbative level. This can be done 
either using the unimproved effective potential in Eqs. (115|17p . or the improved one in 
Eq. (I104p . This latter yields much stronger constraints on the parameter space of the 
theory, thus providing a much finer inspection on the dynamical breaking of the conformal 
symmetry in this model. This fact becomes manifest if we plot sections of the parameter 
space highlighting the region where a valid y could be found. Plots for e 2 = 0.3, 0.6, and 0.9 
are shown in Fig. [5] for the same range of the parameters, the unimproved effective potential 
would pose no restrictions. As an example, for e 2 = 0.9 and n — 1, from Fig. [5] we obtain 
the restriction a > 1.15, so in principle a lower bound 7.99975/x 2 for the mass of the scalar 
is predicted. No such prediction could be made, in this case, using the unimproved effective 
potential. For larger n, this effect is still more dramatic: in Figs.E] and[7] we plot several 
sections of the parameter space, considering the unimproved and the improved effective 
potentials, respectively. 

Another interesting fact is that, for certain values of x, z, and n, Eq. fl 107j) provides 
two viable solutions for y. This is true both for the unimproved as well as for the improved 
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a a a 

Figure 6: Sections of the parameter space of constant e 2 or a, showing where the dynamical 
symmetry breaking occurs, using the unimproved 2-loop calculation of the effective potential. 

effective potential. For example, for e 2 = 0.5, a = 0.5 and n = 1, the unimproved potential 
leads to the equation 

- 200.852?/ 2 + 60.376y - 0.0120817 = , (108) 

for the determination of y, from which we obtain two solutions 

y 1 = 0.30039 , (109a) 
?/2 = 0.00020 . (109b) 

The corresponding masses predicted for the scalar are mi = 7.7907/i 2 and m 2 = 0.00519/i 2 . 
For the same value values of the parameters e 2 , a and n, the improved effective potential 
yields 

-4.75607 x lOV - 4.75649 x 10 8 y 6 - 2.46246 x 10V - 882982.?/ 4 

-24137.9y 3 - 471.335?/ 2 + 60.0824|/ - 0.0379559 = , (110) 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

? 7 ? 

e z e z e z 

e 2 = 0.3 e 2 = 0.6 e 2 = 0.9 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



a a a 

Figure 7: Same as Fig.[6l but using the improved effective potential. It is apparent that when 
e 2 = a, the effective potential is stable for higher values of n; this feature can also be seen in Fig.0 

whose positive and real solutions are 

yi = 0.02540, (111a) 
1/2 = 0.00063, (111b) 

providing mi = 0.18595/i 2 and mi = 0.015843/i 2 . 

Figure[H]depicts the region of the a— e 2 plane, for n = 5, where such a duplicity of solutions 
occurs, both for the unimproved and improved effective potentials. The most important 
difference between the two cases is that the improved effective potential drastically reduces 
the range of parameters where the duplicity happens. Figured shows how the situation 
changes for different values of n, for the second case. 

The pattern in Eqs. (jllO|llip is quite typical: the solution y 2 is smaller than yi. Fixing 
the parameters e 2 = 0.5 and n — 1, y 2 becomes smaller as a increases. At some point, the 
solution yi approaches zero and becomes negative, so it is not counted anymore as a viable 
solution. This behavior is clearly visible at the first graph in Fig. [TUJ For fixed a and n, 
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n = 5 n = 5 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



e 2 e 2 

Figure 8: Regions of the e 2 -a plane, for n = 1, painted according to the number of viable solutions 
y for Eq. (|107p for the unimproved effective potential (left) and for the improved one (right). 
Black, gray and white means two, one, and none solutions, respectively. 

n = l n = 5 n = 10 




Figure 9: Same as Fig. [8] (right), but for different values of n. For larger n, the region where we 
found a unique solution for the conformal symmetry breaking becomes smaller in absolute terms, 
and also in comparison to the region where we found two solutions. 

the situation is reversed: 1/2 becomes smaller as a decreases, as also seen in Fig. [TUJ 

In summary, there are regions of the parameter space of the theory where there are two 
possible vacua, in which the conformal symmetry was broken by radiative corrections. The 
scalar selfcoupling and mass are clearly different for these two vacua. Our numerical studies 
show, however, that for the improved effective potential, the region of the parameter space 
where such a situations takes place is much smaller than for the unimproved potential. 
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ydo -2 ) ydO -2 ) 




Figure 10: Behavior of the two solutions yi and 1/2 (solid and dashed lines, respectively) when 
varying the parameters of the model. 

VII. CONCLUSIONS 

The Renormalization Group Equation is well known to provide better approximations to 
the effective potential of a given model than a pure perturbative calculation up to a given 
loop order. In this work, we pursued the idea of using the RGE to sum infinite subseries 
of the expansion of the effective potential in powers of coupling constants and logarithms 
L = ln(0 2 ///). 

We focused on a Chern-Simons theory coupled to a fermion and a complex scalar field. 
Renormalization group beta-functions and anomalous dimensions should be known up to 
the two-loop order; we collected results already available in the literature and calculated the 
beta-function for the Yukawa coupling and the wavefunction renormalization of the fermionic 
field. With this information, we were able to use the RGE to extract the maximum amount 
of information of the perturbative calculation, obtaining and improved effective potential 
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which, in principle, should allows us to establish more precisely the properties of the model. 
In particular, we were interested in studying the phase where the conformal symmetry 
breaking of the model is broken by the radiative corrections. 

By comparing the outcomes of the standard analysis of dynamical symmetry breaking 
in the model using the standard effective action calculated from loop corrections and the 
improved one, we shown how the latter indeed provides a more precise determination of the 
properties of the model in the broken phase. This should serve as an instructive example 
of the relevance of using the RGE to obtain the maximum amount of information on the 
effective action from a given perturbative calculation. This idea is quite relevant in the 
context of models with classical conformal invariance which is broken at the quantum level, 
for the sake of obtaining the most precise predictions. 

It would be interesting to extend the calculations discussed in this work to higher loop 
orders, to see whether this would imply in some mild refinement of the results presented 
here, or some even more drastic reduction of the parameter space region where the dynamical 
symmetry breaking happens. 
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